hfAgeHfHf2 | R Documentation |
Plots two-stage Hf growth curves in the binary diagram age vs. initial 176Hf/177Hf.
hfAgeHfHf2(which = rownames(Hfinit), RCC = R.crust, xmin = 0, xmax = NULL, ymin = NULL, ymax = DM[,"176Hf/177Hf"], evol.lines1 = FALSE, evol.lines2 = TRUE, violins = FALSE, age.hist = FALSE, kde.bandwidth = 30, rugs = TRUE, pch = labels[which, "Symbol"], col = labels[which,"Colour"], cex = labels[which, "Size"], new = TRUE, main = "Two-stage Hf isotopic growth diagram",...)
which |
character or numeric; names or sequence numbers of the samples to be used |
RCC |
numeric; the 176Lu/177Hf ratio of the intermediate crustal reservoir for calculation of the two-stage model |
xmin |
numeric; minimum of the age axis in Ga |
xmax |
numeric; maximum of the age axis in Ga |
ymin |
numeric; minimum of the 176Hf/177Hf axis |
ymax |
numeric; maximum of the 176Hf/177Hf axis |
evol.lines1 |
logical; plot the evolution lines from present to the |
evol.lines2 |
logical; plot the evolution lines from the |
violins |
logical; overplot the violin plots for each group? |
age.hist |
logical; overplot the frequency histogram of U–Pb ages? |
kde.bandwidth |
numeric; KDE bandwidth. |
rugs |
logical; show the rugs on both axes? |
pch |
plotting symbols |
col |
plotting colours |
cex |
relative size of the plotting symbols |
new |
logical; should be a new plotting window opened? |
main |
main title for the plot |
... |
optional parameters to the underlying function plotWithLimits |
In this case, Hf growth curves in individual samples are plotted solely using a two-stage model (Kemp & Hawkesworth 2003, 2014). Moreover, shown are Hf growth curves for the two main mantle reservoirs, CHUR and Depleted Mantle (DM). The optional small ticks, or rug, on the x axis correspond to two-stage Hf model ages. On the y axis, the extra ticks portray the initial 176Hf/177Hf ratios.
Another option, controlled by the argument age.hist
, is plotting
of the kernel density estimate (KDE) curve for the U–Pb ages along the x axis.
The smoothness of this curve can be controlled by the parameter kde.bandwidth
,
with the 'convenient compromise'
σ = 30
(Vermeesch & Garzanti 2015; Andersen et al. 2018) set as a default value.
In general, the higher is this value, the smoother curve results.
None.
Hf.r
This function is Figaro compatible (except violins). For the Hf growth diagram of age vs. ε(Hf) values see hfAgeEps.
Vojtěch Janoušek, vojtech.janousek@geology.cz
Andersen T, Kristoffersen M, Elburg MA (2018) Visualizing, interpreting and comparing detrital zircon age and Hf isotope data in basin analysis – a graphical approach. Basin Res 30:132–147 doi: 10.1111/bre.12245
Blichert-Toft J, Albaréde F (1997) The Lu–Hf isotope geochemistry of chondrites and the evolution of the mantle–crust system. Earth Planet Sci Lett 148:243-258. doi: 10.1016/S0012-821X(97)00040-X
Bouvier A, Vervoort JD, Patchett PJ (2008) The Lu–Hf and Sm–Nd isotopic composition of CHUR: constraints from unequilibrated chondrites and implications for the bulk composition of terrestrial planets. Earth Planet Sci Lett 273:48-57. doi: 10.1016/j.epsl.2008.06.010
Chauvel C, Garçon M, Bureau S, Besnault A, Jahn BM, Ding ZL (2014) Constraints from loess on the Hf–Nd isotopic composition of the upper continental crust. Earth Planet Sci Lett 388:48-58. doi: 10.1016/j.epsl.2013.11.045
Griffin WL, Wang X, Jackson SE, Pearson NJ, O'Reilly SY, Xu X, Zhou X (2002) Zircon chemistry and magma mixing, SE China: in-situ analysis of Hf isotopes, Tonglu and Pingtan igneous complexes. Lithos 61:237-269. doi: 10.1016/S0024-4937(02)00082-8
Kemp AIS, Hawkesworth CJ (2003) Granitic perspectives on the generation and secular evolution of the continental crust. In: Holland HD, Turekian KK (eds) Treatise on Geochemistry vol. 3, The Crust (ed. R.L. Rudnick). Elsevier-Pergamon, Oxford, pp 349-410. doi: 10.1016/B0-08-043751-6/03027-9
Kemp AIS, Hawkesworth CJ (2014) Growth and differentiation of the continental crust from isotope studies of accessory minerals. In: Holland HD, Turekian KK (eds) Treatise on Geochemistry (Second Edition). Elsevier, Oxford, pp 379-421. doi: 10.1016/B978-0-08-095975-7.00312-0
Lancaster PJ, Storey CD, Hawkesworth CJ, Dhuime B (2011) Understanding the roles of crustal growth and preservation in the detrital zircon record. Earth Planet Sci Lett 305:405-412. doi: 10.1016/j.epsl.2011.03.022
Spencer CJ, Kirkland CL, Roberts NMW, Evans NJ, Liebmann J (2020) Strategies towards robust interpretations of in situ zircon Lu–Hf isotope analyses. Geosci Front 11:843-853. doi: 10.1016/j.gsf.2019.09.004
Vermeesch P, Garzanti E (2015) Making geological sense of 'Big Data' in sedimentary provenance analysis. Chem Geol 409:20–27 doi: 10.1016/j.chemgeo.2015.05.004
hfAgeEps1, hfAgeEps2, hfIso, hfDMage, hfDM2stgAge, hfViolinplot, vioplot
The actual plotting is done by the function plotWithLimits.
sampleDataset("khantaishir_Hf") hfAgeHfHf2(1:39, RCC = 0.022, # Basaltic (Lancaster et al. 2011) evol.lines1 = FALSE, evol.lines2 = TRUE, xmax = 1.5, ymin = 0.2822, col = "blue" , cex = 1.5) hfAgeHfHf2(1:39, cex = 1.5, xmax = 1.3, ymin = 0.2822, evol.lines1 = TRUE, evol.lines2 = FALSE) hfAgeHfHf2(1:39, xmin = 0.35, xmax = 1.3, ymin = 0.2822, ymax = 0.2830, evol.lines1 = FALSE, evol.lines2 = TRUE, age.hist = TRUE, kde.bandwidth = 48, # 'optimal' sigma value (Andersen et al. 2018) cex = 1.5, rugs = TRUE) hfAgeHfHf2(1:39, xmin = 0.35, xmax = 1.3, ymin = 0.2822, cex = 1, evol.lines1 = FALSE, evol.lines2 = TRUE, violins = TRUE)